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Abstract 

We investigate collective synchronous behaviors in random complex networks of limit-cycle os- 
cillators with the non-identical asymmetric coupling scheme, and find a uniform coupling criticality 
of collective synchronization which is independent of complexity of network topologies. Numerically 
simulations on categories of random complex networks have verified this conclusion. 
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1 Introduction 

Traditionally, a network with complex topology is mathematically described by a random graph from 
the ER model proposed by Erdos and Renyi 1 , which has ever dominated scientists' thinking about 
complex networks for almost 40 years. Till recently, the discoveries of small-world networks 2 and scale- 
free networks 3 pioneered a revolution in the theoretic study of complex networks, and, surprisingly, 
many real-life large-scale complex networking systems are in the categories of small- world and(or) scale- 
free networks, ranging from biological to engineering systems, and even to social and economic systems 

mm 

It has been argued that network topologies significantly affect the emergent network dynamics 
and the small-world phenomena and scale-free features of complex networks have led to a fascinating 
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set of common problems concerning how the complexity of network topology facilitates and constrains 
synchronous behaviors of networks jlO|. 

Collective synchronization in a large population of oscillators having natural different frequencies is 
one of the focal points in the synchronization literature, which is a typical phenomenon in the fields of 
biology, physics, and engineering p^j. After the first mathematical description by Wiener [T2| and the 
later fruitful study by Winfree jl3| . Kuramoto refined this connection between nonlinear dynamics and 
statical physics, and formalized the solution to a network of globally coupled limit-cycle oscillators jl4| . 
answering the situation why the oscillators are completely de-synchronized until the coupling strength 
overcomes a criticality Csyn- 

The recent discoveries of topological complexity accelerate the understanding of this cooperative phe- 
nomenon in the frame of random complex networks |15|-|17|. Especially, Hong et al. reported their 
synchronization observations of a larger criticality of coupling strength on small-world networks than 
that of globally coupled networks ^^l- And, the latest investigation stated the absence of critical cou- 
pling strength in frequency synchronization of a swarm of oscillators connected as a scale-free network 
having a power-law exponent 2 < 7 < 3 

All these fruits indicate that the topology of complex networks does play an important role in collective 
synchronization, and different categories of networks may hold significantly different critical coupling 
strengths for the occurrence of synchrony in mutually coupled oscillators. However, it should be pointed 
out that most of these issues hold an assumption that every pair of connected nodes are coupled together 
with the identical coupling strength. In practice, it is more general that pairs of nodes are connected 
with non-identical and asymmetric couplings. Therefore, it is natural to ask the question that whether 
the collective synchronous behaviors of complex networks with non-identical and asymmetric coupling 
schemes still depend on the network topologies? In this paper, we try to explore this question, and the 
answer is that there does exist a uniform criticality of coupling strength to synchronize diversely random 
networks of limit-cycle oscillators. 

The whole paper is organized as follows. Section II describes the phase evolving model of a complex 
network of limit-cycle oscillators, whose couplings are non-identical and asymmetric as specified in this 
work. A uniform critical coupling strength for collective synchrony among oscillators is analytically 
arrived at the same section, whose validity is verified for categories of random complex networks by 
numerical simulations in Section HI. Section IV finally concludes this investigation. 

2 Model description and main result 

We consider a network of N coupled limit-cycle oscillators whose phases 9i,i = 1, 2, • ■ • , N, evolve as 



where Cy is the coupling strength between node (oscillator) i and node (oscillator) j, and is l(or 
0) if node i is connected (or disconnected) with node j. Frequencies uji,i = 1,2,- ■• ,N, are randomly 
distributed following the given frequency distribution g{uj), which is assumed that g{Lu) — g^—Lo). 




(1) 
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C^,^C, = -, z,j = I,--- ,7V, (2) 



Define the non-identical asymmetric coupling scheme 

C_ 

ki 

where C is a positive constant, and ki is the degree of node i, which fits the given degree distribution 
P{k) of a network. Therefore, we have 

N 

dt ' ki 



do j C ^ A ' f n n \ / o \ 



If for every node i, its degree ki = N,i = 1, 2, • • • , N, model © equals the classic Kuramoto model for 
globally coupled networks in |14l IHl E] ■ 

Define the order parameter (r, as 

i^ _ JdLoJdkJ d9g{Lu)P{k)kp{k,Lu;t,0)e'^ 

~ fdkP{k)k ' ^ ' 

where p{k,Lu;t,9) is the density of oscillators with phase 9 at time t with the given frequency uj and 

degree k, which satisfies the normalization as 

p{k,uj;t,9)d9 = 1. (5) 

Assume v to be the continuum limit of the right-hand side (r.h.s.) of Eq. and each randomly 
selected edge is connected to the oscillator having degree k, frequency lu, and phase 9 with probability 
kP{k)g{u)^{k^u],t.e) ^ fpj-^gpgfoi-e, determined by the continuity equation 

dl _ d{pv) 
dt 89 ' 

the density p{k,oj;t,9) evolves as 
dp{k,uj;t,9) d 



dt 89 
Substituting Eq. igj into Eq. (TJl yields 



, ^ C Jdu' J dk' J d9'g{uj')P{k')k'p{k\ Lu'; t, 9') . 



(7) 



8t 89 
whose solution independent of time is 
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^^^'"'"''^'^^ ^~^Ap{k,Lu;t,9)[u + Cr,iniM'~9)]}, (8) 



{p{k, Lu; 9) [uj + Crsin(* -9)]}^ 0. (9) 
where p{k, uj; 9) could be assumed to be 

f 5(6i-arcsin(7^)) if |^ < 1 

^(^'-;^)- ^"^^^ ^J-. (10) 

1 — \ otherwise. 

Here, we assume ■0 = without loss of generality, and D{k,uj) is the normalized factor. Recall g{—Lo) = 
g{Lu), and substitute Eq. pU|) into Eq. I^J. We therefore arrive at 

_ /cifc/_%ci..gMfcP(fc)e'™(^) 

JdkkP{k) ^ ' 

Cr 



Cr 



(12) 



= Cr / cos^9g{Crsm9)d9. (13) 
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From this self-consistency equation we finally come to 

which is surprisingly the same as the critical coupling strength of globally coupled networks investigated 
by Kuramoto |Ol[TT] . 

The main point of criticality H14|l is that, with the non-identical and asymmetric coupling scheme ©J 
different networks may exhibit the same transition of collective synchronous behaviors regardless of the 
significant difference of their degree distributions, i.e., the complexity of network topology. We will verify 
this point through extensive simulations. 



3 Numerical simulations 

Define the average order parameter as 



(15) 



where (• • • ) and [• • • ] denote the averages over different realizations of intrinsic networks and over different 
realizations of intrinsic frequencies, respectively. In all the simulations, rav is averaged over 10 groups of 
networks satisfying the same degree distribution P{k), and each network has 5 sets of frequencies with 
the distribution g{uj). We further specify 

, . f 0.5, if - 1< < 1 
g{Lo) - <^ (16) 

I 0, otherwise. 

Therefore, Csyn « 1-273 by Eq.^ and Eq. (P|l . 

We first numerically simulate the case of scale-free networks of the BA model ^S] whose degree distri- 
bution is in the power-law form of P{k) oc k^^. Figures Ill2l show the simulating outcomes on scale-free 
networks of BA model with m = mg = 3, and N ~ 256, 512, 1024, 2048, respectively. It can be clearly ob- 
served that even having the finite-size effect for N — 256, 512, 1024, 2048, respectively, the average order 
parameter r^^ shows a rapid increase after the critical coupling strength Cgyn = ^^qj^ , indicating that 
when Csyn < Cj the scale-free network of oscillators begin to cluster with synchrony. In Fig. [3 ravN^'"^^ 
shows a similar dependence on the coupling strength C, which rapidly increased when Csyn < C . 

From this simulation example of the same scale-invariant power-law degree distribution P{k) oc k~^ 
with different network scale N , we clearly conclude that the nonzero constant criticality Csyn is in- 
dependent of the finite size of a network, showing a difference from the investigated case of scale-free 
identical-coupled oscillators in |16| . 

To investigate the independence of Csyn on different degree distributions P(fc), we fix the scale N = 
2048, and the average degree (fc) = 6. Therefore, all networks have the same number of nodes and edges 
with different connectivity patterns in the next simulation studies. 

Two main categories of degree distributions of random complex networks are in the forms of Ppower{k) oc 
and Pexp{k) oc e^''. Therefore, we select the famous ER model 1 to generate networks having an 
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Figure 1: The average order parameter rav vs the couphng strength C for scale- free networks of the BA 
model with N ~ 256, 512, 1024, 2048, respectively. All networks are started from m = toq = 3. 




Figure 2: The average VavN^''^^ vs the coupling strength C for scale-free networks of the BA model with 
N — 256, 512, 1024, 2048, respectively. All networks are started from m = mg = 3. 
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Figure 3: The average order parameter Vav vs the couphng strength C for networks generated by the 
ER model (dashed hnc with circle markers), the Local- World model (dotted line with square markers), 
the BA model (solid line with diamond markers) , and the GKK model having power- law exponent 7 = 6 
(dash-dot line with right triangle markers). All networks have the same scale TV — 2048 and the same 
average degree (fc) = 6. 

exponential degree distribution Pexp{k) (x e~'^, and nominate the BA model and the GKK model PU] 
to generate networks having the scale-invariant power-law degree distribution Ppowerik) oc k with 
7 = 3,6, respectively. The proposed local- world evolving network |19j owns a transition between the 
exponential degree distribution and the power-law degree distribution, therefore it is adopted as the final 
prototype for the coming simulations with the parameters M = 10, m = mo = 3. 

The simulating results of Figs. I3I4I are not out of expectation of criticality (|14() . There is a common 
critical coupling strength of the value Cgyn — these four categories of networked limit-cycle 

oscillators, and both the average order parameter rav and the average VavN'^''^^ increased sharply when 
Csyn < C- In other words, the significant difference among complex network topologies does not show 
effect on the collective synchrony of the non-identically asymmetrically coupled limit-cycle oscillators 
0. However, as concluded from 1151 1171 llfij . there is no such a same critical coupling strength of 
collective synchrony in different random complex networks of identically symmetrically coupled limit- 
cycle oscillators. 
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Figure 4: The average ravN^'"^^ vs the coupHng strength C for networks generated by the ER model 
(dashed Une with circle markers), the Local- World model (dotted line with square markers), the BA model 
(solid line with diamond markers), and the GKK model having power-law exponent 7 = 6 (dash-dot line 
with right triangle markers). All networks have the same scale iV = 2048 and the same average degree 



To summarize, we have investigated the collective synchronous behaviors of complex networks of oscil- 
lators having non-identical asymmetric coupling strengths, which show a uniform synchrony criticality 
regardless of the complexity of networking topologies. One of the direct extensions is that for scale-free 
networks, when synchronization is not preferred in some practical situations 21 , the proposed asymmet- 
ric coupling scheme ^ could be applied to obtain a nonzero synchrony criticality which is independent of 
the network topologies. More important to note is the investigation in this paper indicates that, with the 
embedment of complicated coupling schemes, the interactions previously observed between the network 
collective behaviors and network complex topologies are still a small tip of an huge iceberg, most of which 
are left for further explorations. 
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